Comprehensive Characterization of fucAO Operon Activation in Escherichia coli.

Wildtype Escherichia coli cells cannot grow on L-1,2-propanediol, as the fucAO operon within the fucose (fuc) regulon is thought to be silent in the absence of L-fucose. Little information is available concerning the transcriptional regulation of this operon. Here, we first confirm that fucAO operon expression is highly inducible by fucose and is primarily attributable to the upstream operon promoter, while the fucO promoter within the 3'-end of fucA is weak and uninducible. Using 5'RACE, we identify the actual transcriptional start site (TSS) of the main fucAO operon promoter, refuting the originally proposed TSS. Several lines of evidence are provided showing that the fucAO locus is within a transcriptionally repressed region on the chromosome. Operon activation is dependent on FucR and Crp but not SrsR. Two Crp-cAMP binding sites previously found in the regulatory region are validated, where the upstream site plays a more critical role than the downstream site in operon activation. Furthermore, two FucR binding sites are identified, where the downstream site near the first Crp site is more important than the upstream site. Operon transcription relies on Crp-cAMP to a greater degree than on FucR. Our data strongly suggest that FucR mainly functions to facilitate the binding of Crp to its upstream site, which in turn activates the fucAO promoter by efficiently recruiting RNA polymerase.

Research advances in the identification of regulatory mechanisms of surfactin production by Bacillus: a review.

Surfactin is a cyclic hexalipopeptide compound, nonribosomal synthesized by representatives of the Bacillus subtilis species complex which includes B. subtilis group and its closely related species, such as B. subtilis subsp subtilis, B. subtilis subsp spizizenii, B. subtilis subsp inaquosorum, B. atrophaeus, B. amyloliquefaciens, B. velezensis (Steinke mSystems 6: e00057, 2021) It functions as a biosurfactant and signaling molecule and has antibacterial, antiviral, antitumor, and plant disease resistance properties. The Bacillus lipopeptides play an important role in agriculture, oil recovery, cosmetics, food processing and pharmaceuticals, but the natural yield of surfactin synthesized by Bacillus is low. This paper reviews the regulatory pathways and mechanisms that affect surfactin synthesis and release, highlighting the regulatory genes involved in the transcription of the srfAA-AD operon. The several ways to enhance surfactin production, such as governing expression of the genes involved in synthesis and regulation of surfactin synthesis and transport, removal of competitive pathways, optimization of media, and fermentation conditions were commented. This review will provide a theoretical platform for the systematic genetic modification of high-yielding strains of surfactin.

Genomic language model predicts protein co-regulation and function.

Deciphering the relationship between a gene and its genomic context is fundamental to understanding and engineering biological systems. Machine learning has shown promise in learning latent relationships underlying the sequence-structure-function paradigm from massive protein sequence datasets. However, to date, limited attempts have been made in extending this continuum to include higher order genomic context information. Evolutionary processes dictate the specificity of genomic contexts in which a gene is found across phylogenetic distances, and these emergent genomic patterns can be leveraged to uncover functional relationships between gene products. Here, we train a genomic language model (gLM) on millions of metagenomic scaffolds to learn the latent functional and regulatory relationships between genes. gLM learns contextualized protein embeddings that capture the genomic context as well as the protein sequence itself, and encode biologically meaningful and functionally relevant information (e.g. enzymatic function, taxonomy). Our analysis of the attention patterns demonstrates that gLM is learning co-regulated functional modules (i.e. operons). Our findings illustrate that gLM's unsupervised deep learning of the metagenomic corpus is an effective and promising approach to encode functional semantics and regulatory syntax of genes in their genomic contexts and uncover complex relationships between genes in a genomic region.

Genetically encoded transcriptional plasticity underlies stress adaptation in Mycobacterium tuberculosis.

Transcriptional regulation is a critical adaptive mechanism that allows bacteria to respond to changing environments, yet the concept of transcriptional plasticity (TP) - the variability of gene expression in response to environmental changes - remains largely unexplored. In this study, we investigate the genome-wide TP profiles of Mycobacterium tuberculosis (Mtb) genes by analyzing 894 RNA sequencing samples derived from 73 different environmental conditions. Our data reveal that Mtb genes exhibit significant TP variation that correlates with gene function and gene essentiality. We also find that critical genetic features, such as gene length, GC content, and operon size independently impose constraints on TP, beyond trans-regulation. By extending our analysis to include two other Mycobacterium species -- M. smegmatis and M. abscessus -- we demonstrate a striking conservation of the TP landscape. This study provides a comprehensive understanding of the TP exhibited by mycobacteria genes, shedding light on this significant, yet understudied, genetic feature encoded in bacterial genomes.

Reciprocal regulation of enterococcal cephalosporin resistance by products of the autoregulated yvcJ-glmR-yvcL operon enhances fitness during cephalosporin exposure.

Enterococci are commensal members of the gastrointestinal tract and also major nosocomial pathogens. They possess both intrinsic and acquired resistance to many antibiotics, including intrinsic resistance to cephalosporins that target bacterial cell wall synthesis. These antimicrobial resistance traits make enterococcal infections challenging to treat. Moreover, prior therapy with antibiotics, including broad-spectrum cephalosporins, promotes enterococcal proliferation in the gut, resulting in dissemination to other sites of the body and subsequent infection. As a result, a better understanding of mechanisms of cephalosporin resistance is needed to enable development of new therapies to treat or prevent enterococcal infections. We previously reported that flow of metabolites through the peptidoglycan biosynthesis pathway is one determinant of enterococcal cephalosporin resistance. One factor that has been implicated in regulating flow of metabolites into cell wall biosynthesis pathways of other Gram-positive bacteria is GlmR. In enterococci, GlmR is encoded as the middle gene of a predicted 3-gene operon along with YvcJ and YvcL, whose functions are poorly understood. Here we use genetics and biochemistry to investigate the function of the enterococcal yvcJ-glmR-yvcL gene cluster. Our results reveal that YvcL is a DNA-binding protein that regulates expression of the yvcJ-glmR-yvcL operon in response to cell wall stress. YvcJ and GlmR bind UDP-GlcNAc and reciprocally regulate cephalosporin resistance in E. faecalis, and binding of UDP-GlcNAc by YvcJ appears essential for its activity. Reciprocal regulation by YvcJ/GlmR is essential for fitness during exposure to cephalosporin stress. Additionally, our results indicate that enterococcal GlmR likely acts by a different mechanism than the previously studied GlmR of Bacillus subtilis, suggesting that the YvcJ/GlmR regulatory module has evolved unique targets in different species of bacteria.

Multi-Omics integration can be used to rescue metabolic information for some of the dark region of the Pseudomonas putida proteome.

In every omics experiment, genes or their products are identified for which even state of the art tools are unable to assign a function. In the biotechnology chassis organism Pseudomonas putida, these proteins of unknown function make up 14% of the proteome. This missing information can bias analyses since these proteins can carry out functions which impact the engineering of organisms. As a consequence of predicting protein function across all organisms, function prediction tools generally fail to use all of the types of data available for any specific organism, including protein and transcript expression information. Additionally, the release of Alphafold predictions for all Uniprot proteins provides a novel opportunity for leveraging structural information. We constructed a bespoke machine learning model to predict the function of recalcitrant proteins of unknown function in Pseudomonas putida based on these sources of data, which annotated 1079 terms to 213 proteins. Among the predicted functions supplied by the model, we found evidence for a significant overrepresentation of nitrogen metabolism and macromolecule processing proteins. These findings were corroborated by manual analyses of selected proteins which identified, among others, a functionally unannotated operon that likely encodes a branch of the shikimate pathway.

Obtaining Detailed Phage Transcriptomes Using ONT-Cappable-Seq.

Detailed transcription maps of bacteriophages are not usually explored, limiting our understanding of molecular phage biology and restricting their exploitation and engineering. The ONT-cappable-seq method described here brings phage transcriptomics to the accessible nanopore sequencing platform and provides an affordable and more detailed overview of transcriptional features compared to traditional RNA-seq experiments. With ONT-cappable-seq, primary transcripts are specifically capped, enriched, and prepared for long-read sequencing on the nanopore sequencing platform. This enables end-to-end sequencing of unprocessed transcripts covering both phage and host genome, thus providing insight on their operons. The subsequent analysis pipeline makes it possible to rapidly identify the most important transcriptional features such as transcription start and stop sites. The obtained data can thus provide a comprehensive overview of the transcription by your phage of interest.

Cross-regulation of Aps-promoters in Lacticaseibacillus paracasei by the PsdR response regulator in response to lantibiotics.

The PsdRSAB and ApsRSAB detoxification modules, together with the antimicrobial peptides (AMPs)-resistance determinants Dlt system and MprF protein, play major roles in the response to AMPs in Lacticaseibacillus paracasei BL23. Sensitivity assays with a collection of mutants showed that the PsdAB ABC transporter and the Dlt system are the main subtilin resistance determinants. Quantification of the transcriptional response to subtilin indicate that this response is exclusively regulated by the two paralogous systems PsdRSAB and ApsRSAB. Remarkably, a cross-regulation of the derAB, mprF and dlt-operon genes-usually under control of ApsR-by PsdR in response to subtilin was unveiled. The high similarity of the predicted structures of both response regulators (RR), and of the RR-binding sites support this possibility, which we experimentally verified by protein-DNA binding studies. ApsR-P shows a preferential binding in the order PderA > Pdlt > PmprF > PpsdA. However, PsdR-P bound with similar apparent affinity constants to the four promoters. This supports the cross-regulation of derAB, mprF and the dlt-operon by PsdR. The possibility of cross-regulation at the level of RR-promoter interaction allows some regulatory overlap with two RRs controlling the expression of systems involved in maintenance of critical cell membrane functions in response to lantibiotics.

Cross-evaluation of E. coli's operon structures via a whole-cell model suggests alternative cellular benefits for low- versus high-expressing operons.

Many bacteria use operons to coregulate genes, but it remains unclear how operons benefit bacteria. We integrated E. coli's 788 polycistronic operons and 1,231 transcription units into an existing whole-cell model and found inconsistencies between the proposed operon structures and the RNA-seq read counts that the model was parameterized from. We resolved these inconsistencies through iterative, model-guided corrections to both datasets, including the correction of RNA-seq counts of short genes that were misreported as zero by existing alignment algorithms. The resulting model suggested two main modes by which operons benefit bacteria. For 86% of low-expression operons, adding operons increased the co-expression probabilities of their constituent proteins, whereas for 92% of high-expression operons, adding operons resulted in more stable expression ratios between the proteins. These simulations underscored the need for further experimental work on how operons reduce noise and synchronize both the expression timing and the quantity of constituent genes. A record of this paper's transparent peer review process is included in the supplemental information.

The Bacillus subtilis cell envelope stress-inducible ytpAB operon modulates membrane properties and contributes to bacitracin resistance.

Antibiotics that inhibit peptidoglycan synthesis trigger the activation of both specific and general protective responses. σM responds to diverse antibiotics that inhibit cell wall synthesis. Here, we demonstrate that cell wall-inhibiting drugs, such as bacitracin and cefuroxime, induce the σM-dependent ytpAB operon. YtpA is a predicted hydrolase previously proposed to generate the putative lysophospholipid antibiotic bacilysocin (lysophosphatidylglycerol), and YtpB is the branchpoint enzyme for the synthesis of membrane-localized C35 terpenoids. Using targeted lipidomics, we reveal that YtpA is not required for the production of lysophosphatidylglycerol. Nevertheless, ytpA was critical for growth in a mutant strain defective for homeoviscous adaptation due to a lack of genes for the synthesis of branched chain fatty acids and the Des phospholipid desaturase. Consistently, overexpression of ytpA increased membrane fluidity as monitored by fluorescence anisotropy. The ytpA gene contributes to bacitracin resistance in mutants additionally lacking the bceAB or bcrC genes, which directly mediate bacitracin resistance. These epistatic interactions support a model in which σM-dependent induction of the ytpAB operon helps cells tolerate bacitracin stress, either by facilitating the flipping of the undecaprenyl phosphate carrier lipid or by impacting the assembly or function of membrane-associated complexes involved in cell wall homeostasis.IMPORTANCEPeptidoglycan synthesis inhibitors include some of our most important antibiotics. In Bacillus subtilis, peptidoglycan synthesis inhibitors induce the σM regulon, which is critical for intrinsic antibiotic resistance. The σM-dependent ytpAB operon encodes a predicted hydrolase (YtpA) and the enzyme that initiates the synthesis of C35 terpenoids (YtpB). Our results suggest that YtpA is critical in cells defective in homeoviscous adaptation. Furthermore, we find that YtpA functions cooperatively with the BceAB and BcrC proteins in conferring intrinsic resistance to bacitracin, a peptide antibiotic that binds tightly to the undecaprenyl-pyrophosphate lipid carrier that sustains peptidoglycan synthesis.

Nutrient limitation and oxidative stress induce the promoter of acetate operon in Salmonella Typhimurium.

Glyoxylate shunt is an important pathway for microorganisms to survive under multiple stresses. One of its enzymes, malate synthase (encoded by aceB gene), has been widely speculated for its contribution to both the pathogenesis and virulence of various microorganisms. We have previously demonstrated that malate synthase (MS) is required for the growth of Salmonella Typhimurium (S. Typhimurium) under carbon starvation and survival under oxidative stress conditions. The aceB gene is encoded by the acetate operon in S. Typhimurium. We attempted to study the activity of acetate promoter under both the starvation and oxidative stress conditions in a heterologous system. The lac promoter of the pUC19 plasmid was substituted with the putative promoter sequence of the acetate operon of S. Typhimurium upstream to the lacZ gene and transformed the vector construct into E. coli NEBα cells. The transformed cells were subjected to the stress conditions mentioned above. We observed a fourfold increase in the ß-galactosidase activity in these cells resulting from the upregulation of the lacZ gene in the stationary phase of cell growth (nutrient deprived) as compared to the mid-log phase. Following exposure of stationary phase cells to hypochlorite-induced oxidative stress, we further observed a 1.6-fold increase in ß galactosidase activity. These data suggest the induction of promoter activity of the acetate operon under carbon starvation and oxidative stress conditions. Thus, these observations corroborate our previous findings regarding the upregulation of aceB expression under stressful environments.

Exploration of the Native Sucrose Operon Enables the Development of an Inducible T7 Expression System in Paenibacillus polymyxa.

The use of Paenibacillus polymyxa as an industrial producer is limited by the lack of suitable synthetic biology tools. In this study, we identified a native sucrose operon in P. polymyxa. Its structural and functional relationship analysis revealed the presence of multiple regulatory elements, including four ScrR-binding sites and a catabolite-responsive element (CRE). In P. polymyxa, we established a cascade T7 expression system involving an integrated T7 RNA polymerase (T7P) regulated by the sucrose operon and a T7 promoter. It enables controllable gene expression by sucrose and regulatory elements, and a 5-fold increase in expression efficiency compared with the original sucrose operon was achieved. Further deletion of SacB in P. polymyxa resulted in a 38.95% increase in the level of thermophilic lipase (TrLip) production using the cascade T7 induction system. The results highlight the effectiveness of sucrose regulation as a novel synthetic biology tool, which facilitates exploring gene circuits and enables their dynamic regulation.

A type III polyketide synthase cluster in the phylum Planctomycetota is involved in alkylresorcinol biosynthesis.

Members of the bacterial phylum Planctomycetota have recently emerged as promising and for the most part untapped sources of novel bioactive compounds. The characterization of more than 100 novel species in the last decade stimulated recent bioprospection studies that start to unveil the chemical repertoire of the phylum. In this study, we performed systematic bioinformatic analyses based on the genomes of all 131 described members of the current phylum focusing on the identification of type III polyketide synthase (PKS) genes. Type III PKSs are versatile enzymes involved in the biosynthesis of a wide array of structurally diverse natural products with potent biological activities. We identified 96 putative type III PKS genes of which 58 are encoded in an operon with genes encoding a putative oxidoreductase and a methyltransferase. Sequence similarities on protein level and the genetic organization of the operon point towards a functional link to the structurally related hierridins recently discovered in picocyanobacteria. The heterologous expression of planctomycetal type III PKS genes from strains belonging to different families in an engineered Corynebacterium glutamicum strain led to the biosynthesis of pentadecyl- and heptadecylresorcinols. Phenotypic assays performed with the heterologous producer strains and a constructed type III PKS gene deletion mutant suggest that the natural function of the identified compounds differs from that confirmed in other bacterial alkylresorcinol producers. KEY POINTS: • Planctomycetal type III polyketide synthases synthesize long-chain alkylresorcinols. • Phylogenetic analyses suggest an ecological link to picocyanobacterial hierridins. • Engineered C. glutamicum is suitable for an expression of planctomycete-derived genes.

Enhanced Biofilm Formation by Tetracycline in a Staphylococcus aureus Naturally Lacking ica Operon and atl.

Staphylococcus aureus is a major, widespread pathogen, and its biofilm-forming characteristics make it even more difficult to eliminate by biocides. Tetracycline (TCY) is a major broad-spectrum antibiotic, the residues of which can cause deleterious health impacts, and subinhibitory concentrations of TCY have the potential to increase biofilm formation in S. aureus. In this study, we showed how the biofilm formation of S. aureus 123786 is enhanced in the presence of TCY at specific subinhibitory concentrations. S. aureus 123786 used in this study was identified as Staphylococcal Cassette Chromosome mec III, sequence type239 and naturally lacking ica operon and atl gene. Two assays were performed to quantify the formation of S. aureus biofilm. In the crystal violet (CV) assay, the absorbance values of biofilm stained with CV at optical density (OD)540 nm increased after 8 and 16 hr of incubation when the concentration of TCY was 1/2 minimum inhibitory concentration (MIC), whereas at the concentration of 1/16 MIC, the absorbance values increased after 16 and 24 hr of incubation. In tetrazolium salt reduction assay, the absorbance value at OD490 nm of S. aureus 123786 biofilms mixed with 3-(4,5-dimethylthiazol-2-yl)-5-(3-carboxymethoxyphenyl)-2-(4-sulfophenyl)-2H-tetrazolium solution increased after 8 hr when the concentration of TCY was 1/4 MIC, which may be correlated with the higher proliferation and maturation of biofilm. In conclusion, the biofilm formation of S. aureus 123786 could be enhanced in the presence of TCY at specific subinhibitory concentrations.

The alternative sigma factor RpoS regulates Pseudomonas aeruginosa quorum sensing response by repressing the pqsABCDE operon and activating vfr.

Pseudomonas aeruginosa is an important opportunistic pathogen. Several of its virulence-related processes, including the synthesis of pyocyanin (PYO) and biofilm formation, are controlled by quorum sensing (QS). It has been shown that the alternative sigma factor RpoS regulates QS through the reduction of lasR and rhlR transcription (encoding QS regulators). However, paradoxically, the absence of RpoS increases PYO production and biofilm development (that are RhlR dependent) by unknown mechanisms. Here, we show that RpoS represses pqsE transcription, which impacts the stability and activity of RhlR. In the absence of RpoS, rhlR transcript levels are reduced but not the RhlR protein concentration, presumably by its stabilization by PqsE, whose expression is increased. We also report that PYO synthesis and the expression of pqsE and phzA1B1C1D1E1F1G1 operon exhibit the same pattern at different RpoS concentrations, suggesting that the RpoS-dependent PYO production is due to its ability to modify PqsE concentration, which in turn modulates the activation of the phzA1 promoter by RhlR. Finally, we demonstrate that RpoS favors the expression of Vfr, which activates the transcription of lasR and rhlR. Our study contributes to the understanding of how RpoS modulates the QS response in P. aeruginosa, exerting both negative and positive regulation.

BadR directly represses the expression of the glycerol utilization operon in the Lyme disease pathogen.

Glycerol utilization as a carbohydrate source by Borreliella burgdorferi, the Lyme disease spirochete, is critical for its successful colonization and persistence in the tick vector. The expression of the glpFKD (glp) operon, which encodes proteins for glycerol uptake/utilization, must be tightly regulated during the enzootic cycle of B. burgdorferi. Previous studies have established that the second messenger cyclic di-GMP (c-di-GMP) is required for the activation of glp expression, while an alternative sigma factor RpoS acts as a negative regulator for glp expression. In the present study, we report identification of a cis element within the 5´ untranslated region of glp that exerts negative regulation of glp expression. Further genetic screen of known and predicted DNA-binding proteins encoded in the genome of B. burgdorferi uncovered that overexpressing Borrelia host adaptation regulator (BadR), a known global regulator, dramatically reduced glp expression. Similarly, the badR mutant significantly increased glp expression. Subsequent electrophoretic mobility shift assay analyses demonstrated that BadR directly binds to this cis element, thereby repressing glp independent of RpoS-mediated repression. The efficiency of BadR binding was further assessed in the presence of c-di-GMP and various carbohydrates. This finding highlights multi-layered positive and negative regulatory mechanisms employed by B. burgdorferi to synchronize glp expression throughout its enzootic cycle.IMPORTANCEBorreliella burgdorferi, the Lyme disease pathogen, must modulate its gene expression differentially to adapt successfully to its two disparate hosts. Previous studies have demonstrated that the glycerol uptake and utilization operon, glpFKD, plays a crucial role in spirochetal survival within ticks. However, the glpFKD expression must be repressed when B. burgdorferi transitions to the mammalian host. In this study, we identified a specific cis element responsible for the repression of glpFKD. We further pinpointed Borrelia host adaptation regulator as the direct binding protein to this cis element, thereby repressing glpFKD expression. This discovery paves the way for a deeper exploration of how zoonotic pathogens sense distinct hosts and switch their carbon source utilization during transmission.

Insights into the genome of the 'Loco' Concholepas concholepas (Gastropoda: Muricidae) from low-coverage short-read sequencing: genome size, ploidy, transposable elements, nuclear RNA gene operon, mitochondrial genome, and phylogenetic placement in the family Muricidae.

BACKGROUND: The Peruvian 'chanque' or Chilean 'loco' Concholepas concholepas is an economically, ecologically, and culturally important muricid gastropod heavily exploited by artisanal fisheries in the temperate southeastern Pacific Ocean. In this study, we have profited from a set of bioinformatics tools to recover important biological information of C. concholepas from low-coverage short-read NGS datasets. Specifically, we calculated the size of the nuclear genome, ploidy, and estimated transposable elements content using an in silico k-mer approach, we discovered, annotated, and quantified those transposable elements, we assembled and annotated the 45S rDNA RNA operon and mitochondrial genome, and we confirmed the phylogenetic position of C. concholepas within the muricid subfamily Rapaninae based on translated protein coding genes. RESULTS: Using a k-mer approach, the haploid genome size estimated for the predicted diploid genome of C. concholepas varied between 1.83 Gbp (with kmer = 24) and 2.32 Gbp (with kmer = 36). Between half and two thirds of the nuclear genome of C. concholepas was composed of transposable elements. The most common transposable elements were classified as Long Interspersed Nuclear Elements and Short Interspersed Nuclear Elements, which were more abundant than DNA transposons, simple repeats, and Long Terminal Repeats. Less abundant repeat elements included Helitron mobile elements, 45S rRNA DNA, and Satellite DNA, among a few others.The 45S rRNA DNA operon of C. concholepas that encodes for the ssrRNA, 5.8S rRNA, and lsrRNA genes was assembled into a single contig 8,090 bp long. The assembled mitochondrial genome of C. concholepas is 15,449 bp long and encodes 13 protein coding genes, two ribosomal genes, and 22 transfer RNAs. CONCLUSION: The information gained by this study will inform the assembly of a high quality nuclear genome for C. concholepas and will support bioprospecting and biomonitoring using environmental DNA to advance development of conservation and management plans in this overexploited marine snail.

Improving gene set enrichment analysis (GSEA) by using regulation directionality.

To infer the biological meaning from transcriptome data, it is useful to focus on genes that are regulated by the same regulator, i.e., regulons. Unfortunately, current gene set enrichment analysis (GSEA) tools do not consider whether a gene is activated or repressed by a regulator. This distinction is crucial when analyzing regulons since a regulator can work as an activator of certain genes and as a repressor of other genes, yet both sets of genes belong to the same regulon. Therefore, simply averaging expression differences of the genes of such a regulon will not properly reflect the activity of the regulator. What makes it more complicated is the fact that many genes are regulated by different transcription factors, and current transcriptome analysis tools are unable to indicate which regulator is most likely responsible for the observed expression difference of a gene. To address these challenges, we developed the gene set enrichment analysis program GINtool. Additional features of GINtool are novel graphical representations to facilitate the visualization of gene set analyses of transcriptome data, the possibility to include functional categories as gene sets for analysis, and the option to analyze expression differences within operons, which is useful when analyzing prokaryotic transcriptome and also proteome data.IMPORTANCEMeasuring the activity of all genes in cells is a common way to elucidate the function and regulation of genes. These transcriptome analyses produce large amounts of data since genomes contain thousands of genes. The analysis of these large data sets is challenging. Therefore, we developed a new software tool called GINtool that can facilitate the analysis of transcriptome data by using prior knowledge of gene sets controlled by the same regulator, the so-called regulons. An important novelty of GINtool is that it can take into account the directionality of gene regulation in these analyses, i.e., whether a gene is activated or repressed, which is crucial to assess whether a regulon or functional category is affected. GINtool also includes new graphical methods to facilitate the visual inspection of regulation events in transcriptome data sets. These and additional analysis methods included in GINtool make it a powerful software tool to analyze transcriptome data.

Iron limitation indirectly reduces the Escherichia coli torCAD operon expression by a reduction of molybdenum cofactor availability.

The expression of most molybdoenzymes in Escherichia coli has so far been revealed to be regulated by anaerobiosis and requires the presence of iron, based on the necessity of the transcription factor FNR to bind one [4Fe-4S] cluster. One exception is trimethylamine-N-oxide reductase encoded by the torCAD operon, which has been described to be expressed independently from FNR. In contrast to other alternative anaerobic respiratory systems, the expression of the torCAD operon was shown not to be completely repressed by the presence of dioxygen. To date, the basis for the O2-dependent expression of the torCAD operon has been related to the abundance of the transcriptional regulator IscR, which represses the transcription of torS and torT, and is more abundant under aerobic conditions than under anaerobic conditions. In this study, we reinvestigated the regulation of the torCAD operon and its dependence on the presence of iron and identified a novel regulation that depends on the presence of the bis-molybdopterin guanine dinucleotide (bis-MGD) molybdenum cofactor . We confirmed that the torCAD operon is directly regulated by the heme-containing protein TorC and is indirectly regulated by ArcA and by the availability of iron via active FNR and Fur, both regulatory proteins that influence the synthesis of the molybdenum cofactor. Furthermore, we identified a novel regulation mode of torCAD expression that is dependent on cellular levels of bis-MGD and is not used by other bis-MGD-containing enzymes like nitrate reductase.IMPORTANCEIn bacteria, molybdoenzymes are crucial for anaerobic respiration using alternative electron acceptors. FNR is a very important transcription factor that represents the master switch for the expression of target genes in response to anaerobiosis. Only Escherichia coli trimethylamine-N-oxide (TMAO) reductase escapes this regulation by FNR. We identified that the expression of TMAO reductase is regulated by the amount of bis-molybdopterin guanine dinucleotide (bis-MGD) cofactor synthesized by the cell itself, representing a novel regulation pathway for the expression of an operon coding for a molybdoenzyme. Furthermore, TMAO reductase gene expression is indirectly regulated by the presence of iron, which is required for the production of the bis-MGD cofactor in the cell.

Transcription-replication interactions reveal bacterial genome regulation.

Organisms determine the transcription rates of thousands of genes through a few modes of regulation that recur across the genome1. In bacteria, the relationship between the regulatory architecture of a gene and its expression is well understood for individual model gene circuits2,3. However, a broader perspective of these dynamics at the genome scale is lacking, in part because bacterial transcriptomics has hitherto captured only a static snapshot of expression averaged across millions of cells4. As a result, the full diversity of gene expression dynamics and their relation to regulatory architecture remains unknown. Here we present a novel genome-wide classification of regulatory modes based on the transcriptional response of each gene to its own replication, which we term the transcription-replication interaction profile (TRIP). Analysing single-bacterium RNA-sequencing data, we found that the response to the universal perturbation of chromosomal replication integrates biological regulatory factors with biophysical molecular events on the chromosome to reveal the local regulatory context of a gene. Whereas the TRIPs of many genes conform to a gene dosage-dependent pattern, others diverge in distinct ways, and this is shaped by factors such as intra-operon position and repression state. By revealing the underlying mechanistic drivers of gene expression heterogeneity, this work provides a quantitative, biophysical framework for modelling replication-dependent expression dynamics.